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The approach exposed in the recent paper (A. Shebeko, P. Pa- 
pakonstantinou, E. Mavrommatis, Eur. Phys. J. A 27, 143 (2006)) 
has been applied in studying center-of-mass motion effects on the 
nucleon density and momentum distributions in nuclei. We use 
and develop the formalism based upon the Cartesian or boson 
representation, in which the coordinate and momentum operators 
are expressed through the creation and annihilation operators 
for oscillator quanta in the three different space directions. We 
are focused upon effects due to the center-of-mass and short- 
range nucleon correlations embedded in translationally invariant 
ground-state wavefunctions. The latter are constructed in the so- 
called fixed center-of-mass approximation, starting with a Slater 
determinant wave function modified by some correlator {e.g., after 
Jastrow or Villars). It is shown how one can simplify evaluation 
of the corresponding expectation values that determine the dis- 
tributions. The analytic expressions derived here involve the own 
"Tassie-Barker" factors for each distribution. As an illustration, 
numerical calculations have been carried out for the nucleus 4 He 
with the Slater determinant to describe the nucleon (Is) 4 con- 
figuration composed of single-particle orbitals which differ from 
harmonic oscillator ones at small distances. Such orbitals simulate 
somewhat short-range repulsion between nucleons. Special atten- 
tion is paid to a simultaneous shrinking of the center-of-mass 
corrected density and momentum distributions compared to the 
purely (Is) 4 shell nontranslationally invariant ones. 



1. Introduction 

Treatment of the center-of-mass (CM) motion has been 
an attractive subject of exploration in earlier and more 



recent studies of nuclear theory (see, e.g., [1]- [9]). Those 
studies originated from the necessity to remedy a defi- 
ciency of the nuclear many-body wave function (WF), 
namely its lack of translational invariance (TI) wherev- 
er shell-model single-particle (s.p.) WF's are used for 
its construction. This deficiency is important in quite a 
number of cases. 

In the present investigation we adopt the "fixed-CM 
approximation" [10, 11] as a recipe to restore TI of a 
many-body WF which does not possess this proper- 
ty. We apply it when evaluating the elastic form fac- 
tor (FF) F(q) and the nucleon momentum distribution 
(MD) rj(p) for light nuclei, and more specifically for 4 He 
in its ground state (g.s.). Following [9] we prefer to deal 
with the intrinsic quantities which are determined as ex- 
pectation values of appropriate (multiplicative) opera- 
tors that depend on the corresponding Jacobi variables 
and act on the intrinsic WF's. We have seen in [9] that 
the intrinsic density distribution (DD) pi n t(r), being de- 
fined by the Fourier transform of F(q), does not coincide 
with the diagonal part of the one-body density matrix 
(1DM), which is related in a standard manner to the in- 
trinsic MD. In the context, we also note that the term 
"one-body" used here is somewhat conventional. Let us 
mention that F(q) and rj{p) can be related to the differ- 
ent quantities measured via electron-nucleus collisions, 
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respectively, the elastic electron scattering cross sections 
and the inclusive electron scattering cross sections. First 
of all, we mean comparatively simple relations in the 
Born approximation with the plane electron waves. In 
addition, to the so-called approximation of small inter- 
action times (see [12]- [14] and refs. therein) the double 
differential (e, e') reaction cross section becomes propor- 
tional to an integral of rj(j>) over the momentum range 
that is fixed with certain combination (the so-called y 
- scaling variable) of the momentum transfer q and the 
energy transfer uu (cf. [15]). Of course, in the framework 
of these approximations one neglects off-shell effects in 
the electron scattering on bound nucleons and meson 
exchange currents (MEC) contributions to an effective 
electromagnetic (e.m.) interaction with nuclei. The lat- 
ter should be taken into account (see, e.g., refs. [16,17]) 
when describing the electron scattering on nuclei, espe- 
cially at high momentum transfers (in particular, helping 
to remove certain discrepancy between theory and ex- 
periment in the vicinity of the first minimum |F c fc(q)| at 
q 2 = 10/m~ 2 for 4 He). Therefore, any comparison with 
experimental data omitting such physical inputs has a 
restricted character. Nevertheless, in case of light nuclei 
every approximate evaluation of intrinsic quantities, be- 
ing independent of different constraints originated from 
reaction mechanisms, can be compared with microscopic 
("exact") results. In the respect, our addressing to the 
alpha particle seems to be perfectly explicable. 

The aim of this paper is to show to what extent the 
approach developed in refs. [9,12,18,19] can be useful in 
calculations with more realistic WF's than simple har- 
monic oscillator model (HOM) ones. In this connection, 
we consider the CM correction of F(q) and r](p) treated 
on an equal physical footing, viz., by using one and the 
same translationally invariant g.s. WF that incorporates 
the nucleon-nucleon short-range correlations (SRC). One 
should note that despite much interest over the last two 
decades concerning the MD in nuclei [20]- [26], its CM 
correction does not appear to have been properly treat- 



ed except in certain studies, where harmonic oscillator 
(HO) wave functions were used (see, e.g., [12]). Note al- 
so calculations beyond HOM in [27]. The underlying for- 
malism with basic definitions is exposed in the following 
section. Section 3 contains analytic results of our deriva- 
tions beyond HOM, while the corresponding numerical 
results are discussed and compared with experimental 
data in section 4. 

2. The intrinsic form factor, density and 

momentum distributions with short-range 
correlations included 

By definition, the intrinsic (elastic) FF of a nonrelativis- 
tic system with the mass number A and the total angular 
momentum equal to zero is 

F(q) = F mt {q) = ($ mt | exp[tf - (fi - R)] | $ mt ), (1) 

where $i„t is the intrinsic WF of the system (nucleus), 
fi the coordinate operator for nucleon number 1, and 
R = A- 1 Y,t=i h the CM operator. 

In the fixed-CM approximation, according to the 
Ernst, Shakin and Thaler (EST) prescription [11] the nu- 
clear many-body WF with the total momentum P can 
be written in the form: 

| *p) = \P) | <&f n f>, (2) 

where a round bracket is used to represent a vector in 
the space of the CM coordinate, so that \P) means the 
eigenstate of total momentum operator P. The intrinsic 
WF after EST 

(R = | $) 



EST 
int 



) = 



(3) 



[<* \R = 0)(R = 0\ $)]V2 
is constructed from an arbitrary (in general, translation- 
ally non-invariant) WF by requiring that the CM co- 
ordinate R be equal to zero. The corresponding FF is 
the ratio 

A{q) 



F E ST(q) 



A(oy 



A(q) = ($ | (27r) 3 <5(i?)exp[^- (fi - R)] | $). 



(4) 
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Using the Cartesian representation in which 

r = — = (a 1 + a), p = i— = (a 1 — a) 
V2 V2 



(5) 



with the Bose commutation rules, 

[oi,Oj] = 0, [oi,a]] = 5y (i,j = 1,2,3) (6) 

and arbitrary real c-numbers ro and po that meet the 
condition 



roPo = 1, 



one can show (see [9], [19] and Appendix A) that 
A(q)=exp(-^f\u(q), 



U(q) = J dXexpf-^jFiv,^, 



(7) 



(8) 



(9) 



F(v,Z) = ($\0 1 (v + s)0 2 (v)...0 A (if)\$), (10) 
where 

O y (x) — exp(— x*af) exp(fa 7 ) = E^(—x)E 7 (x) (11) 

(7 = 1,..., A) 
with 

*=v^S ff=t-^(X-$) (12) 

and the renormalized "length "parameter 



A-l 



A 



-ro- 



Further, starting from the definition of the intrinsic 
MD (see [9]), 



V(P) = (*<nt I -P/A-p)\ $ int ), 



(13) 



we consider the distribution in the fixed-CM approxima- 
tion: 

($|(2 7 r) 3 ^)^ 1 -J i /A-p)|$) 

Vest(P) = — 5 (14) 

($ | (2 7 r) 3 <5(i?) | $> 



and the Fourier transform 

Vest(p) = (2tt)- 3 J exp(-ipx)N(x)/N(0)dx (15) 
with 

N(x) = ($ | (2vr) 3 u(i?)exp[z(pi - P/A)x\ | $). (16) 

We see the certain resemblance between the structure 
functions N(x) and A(q), viz., both are determined by 
the expectation values of similar multiplicative operators 
with one and the same trial WF <&. Owing to this, using 
the same algebraic technique we get 



N(x) = exp I ) D(x), 



D(x)= I dAexpl-^- 



where 



Po _ „ 
s = =£, v 

and 



irp 
V2A 



F(v',s'), 



(A - iplx) 



(17) 



(18) 



(19) 



Po 



A-l 



A 



-Pa- 



Certain relation of the MD to the corresponding intrinsic 
density matrix has been shown in [9]. 
After this let us assume a trial WF, 



I #) =| Qcorr) = (7(1,2,- •• ,A) | Det) 
with the Slater determinant 



(20) 



Det) 



±= J2 e v V{\ct> Pl {l))---\4> PA {A))}- (21) 



VGSa 



Here e-p is the parity factor for the permutation V, <j) a the 
occupied orbital with the quantum numbers {a} and the 
summation runs over all permutations of the symmetric 
group Sa- 

The ^4-particle operator C — C(f a — rp, p a — pp) 
introduces the SRC and meets all necessary require- 
ments of the translational and Galileo invariance, the 



'Of course, the operator may be spin and isospin dependent 
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permutable and rotational symmetry, etc. However, be- 
ing translationally invariant itself such a model intro- 
duction of correlations does not enable to restore the 
TI violated with such a shell-model WF as the Slater 
determinant. 

What it follows can be used with the Jastrow corre- 
lator [28] 

A 

C=U fiKfi), (22) 

a<P 

where f(r a p) is a two-body correlation factor whose 
deviation from unity occurs only for small distances 
- fp\ less than a correlation radius r c . 
Another popular option goes back to the lectures by 
Villars in [29] (see also [30]) with a unitary operator 



and 



G = cxp(— iG), 



(23) 



G=^g(a,(3), (24) 

a<p 

where the Hermitian operator g(a, (3) acts onto the space 
of the pair (a, (3). In particular, we could follow the sim- 
plest Darmstadt ansatz [31]: 



g(a,P) = -^{s {r a p)p a p + PafjS (fa/3)}, 



(25) 



where s is a function of the relative coordinate r a p = 
r a — rp. Its canonically conjugate momentum p a p — 



\{p*-Pp)- 



Keeping in mind similar constructions we rewrite ex 
pectation (10) as 

F(v,s) = m-v) I Eti-^is) | $(tf)>, 



(26) 



where 

| $(x))=E 1 (x)...E A (x) | $), 

since E^v + s) = E^vjE^s) and [E a (x),Ep( 
(a, P = 1, . . . , A) for any vectors x and y. 
Moreover, we find that 



E(x) r E~\x) = f + 



—= x 

V2 



(27) 



Pa „ 



E(x)pE~ L (x) = p-i-^= x. 



(28) 



Remind that E^ ^ E~ x . In other words, E a (x) is the dis- 
placement operator in the space of nucleon states with 
the label a. 

Due to this property when handling the similarity 
transformation 

C' = E 1 (x)... E A (x)C(f a - r/3, p a - p p )x 

xE^(x)...Ej\x), 
we get 

& = C{E a {x)? a E-\S) - Ep(x)fpE^(x), 

E a (x)$ a E-\x) - E p {x)p>pE-\x)) = 
= C(r a - fp,p a - pp) = C 



i.e. 



C' = C. 



(29) 



Recall that C is a function of all the relative coordinates 
and their canonically conjugate momenta. 
From eqs. (20) and (29) it follows that 

| <S> CO rr(x)} =Ei(x)... E A (x) | $ corr ) = 



= C\ Det{x)). 



(30) 



Here | Det(x)) = E±(x) . . . E A (x) \ Det) is a new Slater 
determinant composed of the renormalized orbitals, 



| <p a (x;a)} = E a (x) \ (j> a {a)) (a = 1,.. ., A), 



(31) 



viz., 



Det{x)) = 



J2 e r V{\cj> pi {x;l))---\4> PA {^A))}. 



V£S A 



(32) 
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In their turn, such orbitals can be evaluated in a concise where according to eq.(A.l) we have the Tassie-Barker 

analytic form as initial ones are linear combinations of Ftb{q) and the HOM FF Fhom(q)- The factor 

the HOM orbitals (see Appendix A). „ , „ + , . „ 

F1] . ' + /s {{l S )±\0E\{^ 2 q)E 1 ^ 2 q)C\{l 3 )±) 

Following (26) we arrive to F^Ja) = — — 

V ' ((Is) 4 | CtC | (Is) 4 ) 

F corr (v, s) = (<S> corr (-v) | Eli-sjE^s) | t> corr (v)) = (38) 

incorporates the SRC in any way. 

/t^ ,/ -\ i At ^t/ ^ /^A i ,/^\ At this point, one can proceed, at least, along the 

= {Det(-v) | E[{-s)Ei{s)C \ Det(v)). (33) 

two guidelines. One of them could be based upon the 
Expressions (8) and (17) with expectations F(v, s) and representation 

F(v ', s'), which are determined by eq. (33), are certain ~ r o - 4 

base for our calculations. ((I s ) I @ -^i( % ^J 7 ^^ 1 ^ 1 ^J 7 ^^ ^ ^ 

2.1 Several working formulae: application to = ((Is) 4 | C\(—q)Ci(q) | (Is) 4 ), (39) 

where 

In special case of the pure HOM (Is) 4 configuration oc- C (~) E ( r ° ~)C(^ ~ )F~ 1 ( r ° ~) 
cupied by the four nucleons in 4 He we have \/2 \/2 

I <f>corr(x)) =| $corr) = C | (Is) 4 ), (34) = ^ ^ + ^2 ^ + | (4Q) 

taking into account that the HOM g.s. | (Is) 4 ) is the Other continuation is prompted by the relation 

vacuum for operators a a (a — 1,...,A). It is the case, , ro ^ - r fr 2 q 2 \ ^ 

i ^ ,-«\ i , i • • i- -i F \(— 1—= q )E 1 (i—= q ) = exp ( — — ) exp(z<f ri), 

where | Det{v)) does not depend on v coinciding with yJ2 V2 \ 4 J 

the initial Slater determinant | (Is) 4 ). Hence, thai gives rise to 

F ls (v, s) = ((Is) 4 | &E{(-3)Ei{3)d | (Is) 4 ). (35) Fcorr{q) = cxp ( r lpj Fc{q)j (41) 



In other words, under such a simplification the function 
F(i 
get 



F(v, s) in integral (9) becomes independent of A and we Fr(ci) — ^—^ — I C exp(z<fri)C | (Is) ) 

W ((Is) 4 I CtC I (Is) 4 ) V ' 

where Fc(q) is the no CM corrected FF with the corre- 



U(q) = U ls (q) j cxp (-^) d\ lated S- S - d I ( ls ) 4 )' 

J \ / A TT ^1 noniicl TITO 

so that 

LUg) _ ((l*) 4 I Ctj^ggA(^gC I (Is) 4 ) 
Uu(0) ~ ((Is) 4 I CtC I (Is) 4 ) 



Analogously, we find 

N(x) = Ntb(x) N H om{x) N corr (x) (43) 



with the own Tassie-Barker factor 

x 

AA 



2 2 x 

(36) 7V TB (x)=exp(^-) (44) 



Thus, the FF of interest is and 
F ES T{q) = F TB (q)F IPM (q)F corr (q), (37) 



N HO M(x)=exp[-^-), (45) 
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((U)*\&E\{fyg)E l {-$ i X)e\{ls) i ) 



V5" 



((la)* | CtC | (Is)*) 



(46) 



Again, different continuations are possible (cf. the 
transition from eq. (38) to eqs. (39) and (41)). In partic- 
ular, with the help of 



^V2^ ^~V2 X ^ = CXP V^T~ ' <,XP 
we get 

N corr (x)=exp(^-)N c (x), 



(tpixj 



N c (x) = 



((Is) 4 | & exp( i # 1 f)C | (Is) 4 ) 



((Is) 4 | CtC | (Is)*) 
The Fourier transform 



Vc(p) 



(2nfj 



e - l P s N c (x)dx 



(47) 



(48) 



(49) 



gives us the one-body momentum distribution (OBMD) 
without the CM correction of the model g.s. C | (Is) 4 ). 
By definition, the intrinsic DD is 

1 



pint(r) 



(2tt) 3 



e-*« r F int (q)dq 



= ($i„t | S(n-R-r) | $ mt ), 
so that the relations 



(50) 



Pest(t) = 



and 



Pc(r) = 



1 



(2tt) 3 



(27T) 3 



e- l9> > BS r(<z)d'f 



e- l ^F c (g)dq 



are, respectively, the one-body density distribution (OB- 
DD) with the CM correction and the no CM corrected 
distribution. 

Thus, we have shown (with the help of purely alge- 
braic means) that evaluation of the distributions can be 



reduced to the well-known treatment. Indeed, expecta- 
tions values (42) and (48) occur in all conventional cal- 
culations with the many-particle WF (20), i.e., without 
any CM correction. Diverse methods have been elaborat- 
ed when evaluating similar quantities (see, e.g., [6], [8], 
[30], [31], [32], [33], [34], [35] and refs. therein). In this 
work we confine ourselves to comparatively simple com- 
putations for a (Is) 4 configuration, where a short-range 
repulsion between nucleons is introduced in an effective 
way, viz., modifying the s.p. orbital as in [10]. Respec- 
tively, the following WF is used in the next section. 

3. Analytic expressions for the form factor, 
density and momentum distributions with 
the single-particle wave function beyond 
HOM 

In accordance with [10] we employ the normalized Rad- 
hakant, Khadkikar and Banerjee (RKB) radial orbital 
for the lowest s.p. state of 4 He, 

1 



<f> RKB (r) 



■■((fao(r) +0<ho(r)), 



(51) 



where <^oo and (fiio are the normalized HO radial eigen- 
functions: 



It* / v" 
<Mr)=2 A / — exp(-— 2 j. 



V^b H b H 



(52) 



<f>w(r) 



3! 



_ 2_f_ 
3^ 



H 



exp 



•55.) < 53 » 



for the states with n — 0, 1 — and n = 1,1 = 0, respec- 
tively. Here bjj is the HO parameter and (3 is a mixing 
parameter. 

The RKB WF allows one to obtain the following ex- 
pressions for the density distribution (normalized to uni- 
ty), for the point proton FF as well as for the MD (also 
normalized to unity): 

p ^ KB{r) = w^hW+W) cxp (" 



6 



A SIMULTANEOUS CENTER-OF-MASS CORRECTION 



(54) 



V6 



24 



(55) 



hi 



H 



tt0F(1 + /3 2 ) 



exp(-6^)x 



-i 2 



(56) 



3.1 The CM corrected form factor F(q) and its 
reduction to quadratures 

Assuming a Slater determinant as the g.s. (Is) 4 of 4 He 
its FF in the fixed-CM approximation can be written in 
the form (cf. [10]): 



F E sT(q) 

where 
Fu(v) = J 



fF ls (\q + u\)F?Ju)du 
fF? s (u)du 



(57) 



/ Avr ±2 



<t> ls (r)dr 



(r) sin vr rdr 



is the no CM-corrected FF. 

In case of the RKB-like s.p. WFs whose orbitals are 
truncated expansions in the radial HO eigenfunctions the 
multiple integrals in the r.h.s. of eq.(63) can be expressed 
through simple integrals. The respective algebraic tech- 
nique has been developed in ref. [12] and exposed recent- 
ly in [37] (see also Appendix A to the present paper). Its 
application with the RKB orbital enables us to get 

A RKB (q) 



? RKB 



A RKB {oy 



(58) 



where 



\RKB 



(q) = h{q) + h(q), 



M 9 ) = ^exp(-A^) 



j{B 2 

o 



Ml 



!/ 1 , si 
- A {t+-qb H f 



\{t+\qb H f 



} enp{-t 2 )tdt, 



h(q) = ^-exp(-^q 2 b 2 H ) > 



o 



Ml 



!/ 1 , S2 



+ 



+B 2 



\{t + \qb H f 



Ml 



\(t-\qb H f 



} exp(-t 2 )dt, 



i RKB 



(0) 



oc 

4tt r 



M 2 4 



if- 



exp(-t 2 )rrft 



The functions M 2 (z) and B 2 (z) are second degree poly- 
nomials of the variable z 



- RKB+CMM 
Frosch et al. 




Fig.l. The point-like FF (left) and the charge FF (right) of A He. 
Curves calculated with RKB WF using the EST prescription (sol- 
id) and without the CM-correction (dashed); experimental points 
from [39]. Other clarifications in the text 



M 2 {z) = mo + m\z + m 2 z 2 
B 2 (z) = h Q + h x z + h 2 z 2 , 

where the constants are related to the mixing parameter 



mo = 1+/3 2 , 

mi = 2^/273/3(1 - y/2/3/3), 
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m 2 = (2/3)/3 2 , 
and 

ho = mo + mi + 2m 2 , 
hi = mi + 2m 2 , 
h 2 = m 2 . 
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Fig. 2. The one-body density distribution (OBDD) for (Is) 4 config- 
uration with RKB orbital in the fixed-CM approximation (solid) 
and without the CM correction (dashed). For two sets of parame- 
ters: b H = 0.8532 fm and = -0.4738 (top); b H = 0.8532 fm and 
/3 = (bottom). Dot-dashed with the parametrization from [41] 

3.2 The CM corrected momentum distribution 
r)(p) and its reduction to quadratures 

In parallel, starting from eq.(16), we obtain for the (Is) 4 
configuration with the Slater determinant | $) =| (Is) 4 ) 
(see Appendix B to Lect.I in [37]): 

Nest(x) — j dk (Is | exp ( j exp (^~~^~P^j I ^ s ) x 



p (fm') 



Fig. 3. The one-body momentum distribution (OBMD) for the 
(Is) 4 configuration with RKB orbital. Difference between the 
curves is the same as in fig. 2. Dot-dashed with parametrization 
from [42] 

Again, using the representation (5) and splitting the 
exponents involved in eq.(67) with the succesive normal 
ordering of the operators a t and a (the former are to 
the left from the latter), one can show that 

4-1 i 2 \ J(x) 



N{x) RKB _ 

N (Q)RKB ~ CX P ^ A 4b 2 H J J(0) • 

where the integral J(x) is determined by 



(60) 



J{x) = j 



exp - 



4A 



g(X 2 ;x 2 )X 2 d\ 



g(X 2 ;x 2 ) = + i.l.) Bf + (.1, + j B'{B 2 + 
+3 QAi + ^A 2 \ B\B\ + (^A x + X -A 2 ) B\ 



A l = 1 + I3 2 



1" 



A-l\ 2 P 2 x 2 



+ 



2A 2 



x (Is | exp I i— I exp 



px 



I l^) 3 . 



(59) 



A-l\ 2 p 2 x 2 



rfA 2 
2A 2 
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2 z 2 A 2 



A 2 



B 1 = l + p 2 -\ -(3 



2 2 

2A 2 



+ 



/3 



2A 2 



vlx 1 _ rp?_ 
2A 2 2A 2 



Bo 



6 P ^l 4 



Thus, the structure function N RKB (x) can be reduced 
to one-dimensional integrals similar to those derived for 
F RKB (q). Here A = 4, but we allow A to be changeable, 
particularly, in order to check that the corresponding 
distribution 



RKB 

Vest 



(p) = — L- J N RKB (x)/N RKB (0)sm(px)xdx 
o 



to the limit A — > oo yields the no CM corrected distri- 
bution (62). 



4. Results and discussion 

Analytic expressions obtained in sect. 2 for the density 
and momentum distributions and their Fourier trans- 
forms are sufficiently general to be applied in different 
translationally invariant treatments with the SRC in- 
cluded. The corresponding formulae derived in sect. 3 in 
case of the 4 He nucleus have been employed to carry out 
our calculations beyond the simple HOM. Their numer- 
ical results are displayed in figs. 1-4. 
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Fig. 4. Variations log rj^^(p) (solid), log rj^ <B (p) (dashed) and 
log r/Moritaip) (dot-dashed) with p. Notation log riMorita(p) = 
(f ) W SAr (fp)> where function W SN (x) calculated with a conve- 
nient parametrization from [42]. Points are resulted from [24] 



In fig. 1 we show our calculation of the charge FF, 
Fch(q 2 ) = fp{ ( l)F{q) of the alpha particle, using eq.(55) 
and eq.(58) and considering for the finite proton size fac- 
tor f p (q) the Chandra and Sauer prescription [40]. The 
two parameters bu and (3 have been determined by the 
least-square fitting to the experimental values [39] : their 
best-fit values are b H = 0.8532 fm and (3 = -0.4738 
(\ 2 ~ 13.07). These values have been utilized in our 
calculations shown in figs. 2-4. 

As it is seen in fig.l (its left part), the CM cor- 
rection leads to a considerable qualitative change of in- 
dependence of the FF: its first minimum and second 
maximum are shifted towards higher q- values. This dif- 
ference between solid and dashed curves is due to the 
different behavior of the respective densities at small dis- 
tances r < 1 fm: see fig. 2, top, where the dashed curve 
is p R ^ B {r) by eq.(54), while the solid is determined by 



RKB 



(r) 



2vr 2 r 



F RKB (q) sm(qr)qdq. 



(61) 
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Moreover, we see that each of the distributions (for the 
simple HOM orbital on the bottom and for the RKB 
orbital on the top), after being CM corrected, increases 
in its central but decreases in its peripheral region. One 
may say that we encounter a specific effect of shrinking 
the OBDD owing to the translationally invariant treat- 
ment. 

In addition, there is a central depression of the den- 
sity distribution (cf. the upper and lower dashed lines in 
fig. 2). Such a change is not unexpected since the RKB 
WF represents a simple way to allow for some of the 
effects of short-range repulsion between the nucleons in 
4 He. These numerical results get an explicit confirma- 
tion if one writes 



RKB 



(62) 



Evidently, the inequality 



HK B ^ HKB( n \ I _ HOM / n \ 

Psp (0) < P sp (0) 1/3=0 = p (0) 



takes place for negative (3 values with | (3 |< 2\/6- 

In parallel, we show in fig. 3 that the corresponding 
change of the OBMD has much in common with that for 
the OBDD, viz., the distribution t]est(p) turns out to 
be shrunk in the above sense relative to the distribution 
7 l^ p KB (p)- Thus, we see a simultaneous shrinking of the 
density distribution p{r) and the momentum distribu- 
tion rj(p). As has been shown in [12] (see also [14]), such 
a simultaneous change of these distributions plays a sub- 
stantial role in getting a fair treatment of the data on 
the elastic and inelastic electron scattering of 4 He. Let 
us recall that there the charge FF and the dynamic FF 
of 4 He were calculated using one and the same HOM 
WF, corrected both with the fixed-CM approximation 
and the Peierls-Yoccoz prescription [44]. 

Regarding the properties of these simultaneously cor- 
rected distributions in detail, we would like to emphasize 
a practical consequence of their interpretation. This as- 
pect becomes especially transparent in the case of the 



simple HOM (Is) 4 configuration, where we have 

P^S 4 (r) = {^f }- 3 cM-r 2 /ft) 



vs 



HOM 

rsp 



(r) 
and 



[V^r ] 3 exp(-r 2 /ro) 



vs 



HOM 
fsp 



(p) = W^Po] 3 cxp(-p 2 /Po) 



Thus, the inclusion of CM corrections gives rise to the 
two independent renormalizations, tq = bn — > fo = 
y/Z/ 4ro and po = bn 1 — > Po = ^/3/4po, of the oscil- 
latory parameter values, tq and po (cf. [12]). Evidently, 
such changes are not equivalent to a hasty replacement 
of po by ^/4/3po if one follows the Tassie-Barker recipe 
with bn — * y/Z/Abn only. 

Now, following the conventional way of determining 
the HOM parameter tq, as in ref. [12], we will use the 
expansions 



F ch (q 2 ) = 1 



1 



2 2 
1 r ch + ' 



2 2 

<rK + ■ 



f P (q) 
and 

F(q)- 

where we have in HOM 



6" p 

1 1 2 2 

6 



2 ° 

<p = f 

rms 2 1 

so that 



r ch — 2 r + T p ' 

whence 



„2 2 



'0 



exp 



(63) 



Doing so for the CM corrected quantities we find the 
similar relation 

o 2 r o ti o 



' ch 



exp 



(64) 
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with the identical g-dependence Fest 1 (l) = 



?HOMi 



(q) = exp(— q r e x „/4). At the same time the dif- 
ference between the respective OBMDs becomes more 



considerable than after the substitution r — > <J -j^r exv 
in r]" P OM (p) = exp(-p 2 rg) that gives 



HOM 
'Isp 



x exp 



vs 



(p) = 



A 



A-l 



2 2 



3/2 r 3 



r 3/2 



sp 



(p) 



T 3/2 



exp(-p 2 r 2 ). 



Under the simultaneous CM correction of the OBDD and 
OBMD we have 

A 



A-l 



3 „3 
exp 



r 3/2 



x exp 



vs 



HOM i 



A 



A-l 



2 2 



(65) 



(66) 



that is equivalent to the substitution ro — > -j^j 

Note also that the product ropo = 1 — ^4 _1 7^ 1, unlike 
the relation r^po = 1. In this connection, following [19] 
let us remind the commutation rules for intrinsic coordi- 
nate r ' = r — R and conjugate momenta p' = p — P/A 



[r;,p']=^u(l-lM) ; {1,3 = 1,2,3) 



(67) 



One can show that the corresponding uncertainty princi- 
ple is related to the deviation from unity. Thus, the un- 
certainty principle does not contradict the simultaneous 
shrinking of the density and momentum distributions 
(see also [37], Lect.I, Suppl. C) 

In case of the RKB function we get 



2 — _ 2 

^rras 2^ 



JV6_ 

l + (3 2 



1 




(68) 



It means that the short-range repulsion involved in the 
WF with a negative j3 leads to some increasing the rms 
radius, viz., r™ s B > r^ s M . For the values b H = 0.8532 
fm and = -0.4738 the formula (68) yields r^ s B = 
1.429 fm, so that the corresponding charge radius is 
equal to = 1.667 fm. Here we employ the charge 

proton radius r p = 0.86 fm (see, for example, Appendix 
7 in [43]). 

The CM correction gives an opposite effect. Indeed, 
after some calculation we find that r^fj = 1.309 fm for 
the same b H = 0.8532 fm and (i = -0.4738. From this 
it follows that r B ^ T = 1.566 fm. 

The variation of log r]§g B and of log r]^ 3 with p is 
depicted in fig. 4 for a wider range of momenta. It is seen 
from fig. 4 that the allowance of the CM motion improves 
the description of the available data on the OBMD of 
the alpha particle. It is further seen from fig. 4 that in 
the translationally invariant quantity according to the 
fixed-CM prescription the "seagull" behavior appearing 
in the variation of the corresponding s.p. one becomes 
somewhat less pronounced. The dip is diminished and it 
moves to smaller values of momentum. 

Finally, we would like to point out in connection to 
the comparison with the s.p. distributions that the CM 
corrected OBDD and OBMD become closer to the cor- 
responding microscopic ones by using their convenient 
parametrizations from [41] and [42], as one can see in 
figs. 3 and 4. According to the communication [42] one 
has to introduce the factor (2n)~ 3 to reproduce the 
momentum distribution rjMoritaip) which is one of the 
significant results obtained by the Sapporo group. At 
this point, let us recall that these authors employed the 
so-called ATMS-method, where ATMS is abbreviation 
"Amalgamation of Two-body correlations into Multiple 
Scattering process", to construct the variational WF of 
the He nucleus (see [21] and refs. therein). Along the 
variational approach a considerable progress was made 
when including more dynamics of the realistic nucleon- 
nucleon interaction such as the effect of its tensor com- 
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ponent (cf. [22]). 



5. Concluding remarks 



Also, it is a great pleasure for us to thank A. Antonov 
and M. Gaidarov for sending the numerical values of the 
momentum distribution displayed by the points in fig. 4. 



We have seen how the approach exposed in [9] can be 
extended to the translationally invariant evaluation of 
the density and momentum distributions in nuclei. The 
present analysis shows that the restoration of transla- 
tional invariance in the Slater determinant WF of i He 
by means of the fixed-CM correction (the EST prescrip- 
tion) gives rise as a whole to essential changes in the 
r-, p- and ^-dependences of the OBDD p(r), the OBMD 
rj(p) and the charge FF F c h(q), respectively. We have 
seen that the correlation between nucleons induced by 
the fixation of the center-of-mass of the nucleus results 
in the simultaneous shrinking of p(r) and rj(p). Mean- 
while, this effect has been revealed here beyond the pure 
HOM extending the available experience. 

Also, this study demonstrates the relative impor- 
tance of the CM and SRC corrections for the same nu- 
cleus, viz., the shrinking of the density and momentum 
distributions owing to the use of translationally invari- 
ant g.s. wave functions of 4 He vs their broadening af- 
ter the inclusion of short-range repulsion in these wave 
functions at small distances r < 1 fm. It is true that 
the latter has been introduced in our calculations in 
a simple manner. Nevertheless, there are all reasons to 
believe that the algebraic method employed here might 
be helpful within more sophisticated approaches, where 
the short-range correlations are taking into account via 
the Jastrow factor or other correlation operator (see, 
e.g. [31]). At present, the corresponding applications are 
in progress both for the A He and 16 nuclei. 
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Appendix A 

Some details of calculations beyond HOM 

Here we want to illustrate a convenient method for 
evaluation of the expectations in question being aimed at 
some general (model-independent) results (cf. [9], [19]). 

First of all, we have by recurring the Cartesian rep- 
resentation: 



exp 



iq (fi - R) 



exp 



iq 



A-l 



exp 



~ iq A 



= F TB (q) F HO m (q) x 
'A-l\ r 



x exp 



x cxp 



iq 



iq 



A 
A-l 



V2 
V2 



x exp 



-iq 



V2A 



Ol 



exp 



-. r 
-iq 



V2A 



«2 



(Al) 



F T B(q) = exp(^|-), F HOM (q) = exp(-^-), where 
the index a at a a (a|J is the individual particle number 
(a = l,---,A). 

Thereat, the Tassie-Barker factor FrB(q) appears au- 
tomatically due to a specific structure of the operators 
involved. In other words, its appearance is independent 
of any nuclear properties (in general, properties of a fi- 
nite system). The only mathematical tool that has been 
used is the Baker-Hausdorff relation: 



.A+B 



e A e B e -|[A,B] 



(A2) 



that is valid with arbitrary operators A and B for which 



the commutator 



A, B 



commutes with each of them. 



Further, applying eq.(A.2) in combination with 



(2tt) j 5 



(M) = J exp ( 



i\Rj dX, 



(A3) 
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one can show that the expectation value A(q) in eq. (4) 
and the expectation value N(x) in eq. (16) are expressed 
through one and the same function F(x, y) that depends, 
respectively, on the arguments x = v, y = s (as in eq.(9)) 
and x = v', y = s' (as in eq.(18)). In other words, we 
have constructed the common generating function for 
each of them. One should stress that this result has been 
obtained independently of the model WF $. 

The algebraic technique shown here turns out to be 
useful for practical calculations with the Slater deter- 
minants like | <&) (see [18]) or the Slater determinants 
modified by different correlators (for instance, the Jas- 
trow factor). 

In the simplest case of the independent particle mod- 
el (IPM) (Is) 4 configuration for 4 He with the Slater de- 
terminant | $) =| (Is) 4 ) we get omitting the nonessential 
factor [A!]" 1 , 



\IPM 



(q) = exp 



U 



IPM 



(?), 



u 



IPM 



(.)-/ 



dX exp 



AA 



f(X,q), 



/(A, q) = (Is | exp (—ctcft) exp (a a) | Is) x 



x(ls | exp(-(3*tf)exp{/3a) | ls)' c 



(AA) 
(A5) 

(A6) 



[X + (A-l)q], /?: 



i^2-[X-q\, 



V2A V2A 
with the renormalized "length "parameter 



r 



1 



7-0 



In our case ro = bn ■ We are keeping the notations with 
a mass number A (= 4) to point out certain trend in the 
^-dependence. 

At this point, let us note that the RKB orbital (or an- 
other model orbital) being composed of the basis states 
of the spherical representation can be written then as a 
superposition of the basis states | 12^2113} of the Carte- 
sian representation (see, e.g., the monograph [38] and 
refs. therein), 

1 



Vn-i!n 2 !n 3 ! 



(^l t )" 1 (^ t )" 2 (^3 t )" 3 |0>, 



(A.7) 



where the vector | 0) =| 000) is the vacuum state with 
respect to the destruction operators cii (i — 1,2,3), e.g., 

k I 0) = 0. (A8) 

It is proved that for the RKB-orbital, 

I Is) = [1 + f3 2 ]-^ 2 [l - [fi/y/%) a tS t ] 1 0). (A9) 

Substituting (A. 9) into (A. 6) (when calculating the ratio 
A IPM (q)/A IPM (0), the normalization factor [1+/3 2 ]" 1 / 2 
can be omitted) we find 

exp (x • 3) I Is) = [1 - {0/V6)(8 f + x)(a t + ft] \ 0) 

(A10) 

for any complex vector \. 

Now, after modest effort we obtain 

(Is I exp (-x * • a f ) exp (x • a) \ Is) = 
= l + /3 2 -^/? 2 X*X- 

K *x * + x x } + y & *)(* $ ( AJ1 ) 

From (A. 11) it follows, for instance, 

(Is I exp {id* a t) exp (iaa) \ Is) = M 2 . ( A - 12 ) 

where the polynomial M<i(z) is given at the end of sub- 
sect. 3.1. 
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